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Abstract. Two classes of gamma-ray bursts have been identified so far, characterized by durations shorter and 
longer than approximately 2 seconds. In 1998 two independent papers indicated the existence of the third class of 
the bursts roughly duration between 2 and 10 seconds. In this paper, using the full BATSE Catalog, the maximum 
likelihood estimation is presented, which gives a 0.5% probability to having only two subclasses. The Monte-Carlo 
simulation confirms this probability, too. 

Key words. Gamma rays: bursts, theory, observations - Methods: data analysis, observational, statistical 



1. INTRODUCTION 



In the BATSE Current Catalog flMeegan, et al. 200l|) 
there are 2702 Gamma-Ray Bursts (GRBs), of which 2041 
have duration information. Kouveliotou, et al. 1992 have 



identified two types of GRB based on durations, for which 
the value of Tgo (the time during which 90% of the fluence 
is accumulated) is smaller or larger than 2 s, respectively, 
and exhibits an acceptable bimodal log- normal ("two- 
Gaussian") fit. This bimodal distribution has been further 
quantified in another paper ( Kouveliotou, ct al. 1995| ), 
where a two-Gaussian fit was made, however the best pa- 



rameters of the fit were published in [McBreen, et al. 
and 



Koshut, et al. 1996 



1994 



Previously we have published an article 



(Horvath, 1998), where both two- and three-Gaussian fits 
were made using the % 2 method, which gave a 99.98% 
significance the third Gaussian is needed. This is an 
agreement with the result of Mukherjee, et al. 1998, 
who used a multivariate analysis and found that the 
probability of existence of two groups, rather than 
three ones, is less than 10~ 4 . Hakkila, et al. 2000c also 
confirmed this result by statistical clustering analysis. 
However, they suggested that the third group was 
caused 



by 



Hakkila, et al. 2000b) 



instrumental biases 
Recently 



(|Hakkila, et al. 2000a 



alastegui et al. 2001 



have applied automatic classifier algorithms and obtain 
three different classes of GRBs. Add also that the 
intermediate subgroup shows a remarkable angular distri- 



bution on the sky (Balazs et al. 1998, 


Balazs et al. 1999 


Meszaros ct al. 2000 




Litvin ct al. 2001 


)■ 



Although the high probabilities for the occurrence of 
third (intermediate in duration) subgroup are suggestive, 




log T90 



Fig. 1. Duration distribution of the observed BATSE 
bursts. 



the existence of the third subgroup is a matter of debate. 
Hence further studies concerning this subclass are highly 
needed. In order to make further progress in quantify- 
ing this classification, one of the issues which needs to be 
addressed is an evaluation of the probabilities associated 
with the bimodal, or - in general - with the multimodal 
distribution. In this paper we take another attempt at the 
trimodal log- normal distribution of Tgo, evaluating a new 
small probability of the assumption that the third subclass 
is a chance occurrence. 

In Section 2 bi and a tri-modal log-normal fits have 
been made using the maximum likelihood method. In 
Section 3 a hundred Monte-Carlo simulations have been 
taken confirm the low probability. In Section 4 a triggering 



I. Horvath: Further Study of the Gamma-Ray Bursts Duration Distribution 



Table 1. Two Gaussian fit of the GRBs duration distri- Table 2. Three Gaussian fit of the GRBs duration distri- 



bution 



buti 





Duration(logTgo) 


a(logTgo) w 


short 


-0.11 


0.61 0.32 


long 


1.54 


0.43 0.68 



Averagcduratioii(logTQo )standarddcviation(fr) and weight (w)ofthcgroups. 





Duration(logTgo) 


o(logT 90 ) 


w 


short 


-0.25 


0.53 


0.26 


long 


1.55 


0.42 


0.68 


intermediate 


0.63 


0.20 


0.06 



Averagcduration(logTgQ )standarddcviation(cr)andwcight (w)of'thegroups. 



systematic effect has been discussed. Finally the conclu- 
sions are given in Section 5. 

2. FITS IN logT 90 

For this investigation we have used a smaller set of 1929 
burst durations in the Current Catalog, because only they 
have peak flux information as well. We use the Tqq mea- 
sures provided in this data set. Figure 1. shows the distri- 
bution of log Tqq . 

A fit to the duration distribution has been taken using 
a maximum likelihood method with the superposition of 
two log-normal distributions. This can be done by a stan- 
dard search for 5 parameters with N = 1929 measured 



points (cf. Press et al. 1992; Chapt. 15). Both log-normal 
distributions have two parameters; the fifth parameter de- 
fines the weight (w\ ) of the first log- normal distribution. 
The second weight is t«2 = (1 — w\) due to the nor- 
malization. Therefore we obtain the best fit to the 5 pa- 
rametcrs through a ma ximum likelihood estimation (e.g., 
Kendall fc Stuart 1976j ). We search for the maximum of 
the formula 



N 

L = J^ln (w 1 f 1 (x i ,T 1 ,a 1 ) 



W2jl[Xi 1 T- 1 ,(J2)) 



where 



fk 



(Tfc V2tt 



exp 



(i) 



(2) 



where Tj. is the mean in log T go and a is the stan- 
dard deviation. This fit gives us the best parameters of 
the two-Gaussian fit (Table 1.), which are very similar to 
previously published values ( Horvath, 1998[ ). 

Secondly, a three-Gaussian fit has been taken with 
three fk functions with eight parameters (three means, 
three standard deviations and two weights). For the best 
fitted parameters see Table 2. The best logarithm o f 
the likelihood (L 3 ) is 12326.25 ( |Kcndall fc Stuart 1976j ). 



For two Gaussians the maximum of the likelihood was 
L 2 =12320.11. According to the mathematical theory, 
twice the difference of these numbers follows the % 2 dis- 
tribution with three degrees of freedom because the new 
fit has three more parameters 



2(L 3 - L 2 ) 



xt, 



(3) 



The difference is 6.14 which gives us a 0.5% probabil- 
ity. Therefore the three-Gaussian fit is better and there is 
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Fig. 2. Distribution of the MC simulated L3 — L 2 . L3 is a 
likelihood with three-Gaussian and L 2 is a likelihood with 
two-Gaussian. 

a 0.005 probability that it is caused by statistical fluctua- 
tion. 

3. MONTE-CARLO SIMULATION 

One can check this 0.005 probability using the Monte- 
Carlo (MC) simulation. Take the two-Gaussian distribu- 
tion with the best fitted parameters of the observed data, 
and generate 1929 numbers for Tgo whose distribution fol- 
low the two-Gaussian distribution. Find the best likeli- 
hood with five free parameters (two means, two disper- 
sions and two weights; but the sum of the last two ones 
must be 1929). Secondly make a fit with three-Gaussian 
distribution (eight free parameters, three means, disper- 
sions and weights). Take a difference between the two log- 
arithms of the maximum likelihood, which gives one num- 
ber. 

This procedure is repeated 99 times and we have a hun- 
dred MC simulated numbers. Only one of these numbers 
is bigger than the value obtained from the BATSE data 
(6.14). The distribution of these differences are given in 
Figure 2. Therefore the MC simulations confirm the like- 
lihood law statement and gives us a similar probability 
(1%) that the third group is a statistical fluctuation. 

Also the population of the third group generated by 
MC simulations is far from the GRB third group popula- 
tion which is 6% (see Table 2.). The average of the flue- 



I. Horvath: Further Study of the Gamma-Ray Bursts Duration Distribution 



3 



tuation population is 2,5% and only one of the hundred 
number is bigger than 6%. 

4. SYSTEMATICS 

In this section the possibility that the third intermedi- 
ate subgroup is an instrumental effect is discussed. The 
BATSE on-board software tests for the existence of bursts 
by comparing the count rates to the threshold levels for 
three separate time intervals: 64 ms, 256 ms, and 1024 ms. 
The efficiency changes in the region of the middle area 
because the 1024 ms trigger is becoming less sensitive as 
burst durations fall below about one second. This means 
that at the "intermediate" timescale a large systematic 
deviation is possible. To reduce the effects of trigger sys- 
tematics in this region we truncated the dataset to include 
only GRBs that would have triggered BATSE on the 64 
ms timescale. 

Using the Current BATSE catalog CmaxCmin table 
(Meegan, et al. 2001) we choose the GRBs, which num- 
bers larger than one in the second column (64 ms scale 
maximum counts divided by the threshold count rate). 
This process reduced the bursts numbers very much, 
therefore unfortunately just 958 bursts satisfied the above 
condition. 

We repeated the maximum likelihood test with these 
burst's durations. The computed probability is still below 
1%. Therefore after eliminating some systematic effects, 
the third group is still statistically significant. 

5. CONCLUSIONS 

1. It is possible that the three log-normal fit is accidental, 
and that there are only two types of GRBs. However, 
if the Tg distribution of these two types of GRBs is 
log-normal, then the probability that the third group 
of GRBs is an accidental fluctuation is less than 0.5 %. 

2. Therefore, statistically the third component exis- 
tence is not questionable. However, the physical ex- 
istence of the third group is still argued. The sky 
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as proven by 


Balazs et al. 1998, 


Balazs et al. 1999. 


Meszaros et al. 2000 




Litvin et al. 2001 


The logN- 



logS distribution is may also differ from the 



other group distribution Horvath, 1998. Opponently 
Hakkila, et al. 2000c believe the third statistically 
proved subgroup is only a deviation caused by a com- 
plicated instrumental effect, which can reduce some 
faint long burst's duration. This paper does not deal 
with the effect mentioned above, however the trigger- 
ing systematic effects are examined and after that the 
third group is still statistically significant. 
Therefore, this theme should be discussed in future 
papers to further elucidate the reality and properties 
of the third class. 



Acknowledgements. This research was supported in part 
through NATO advanced research fellowship 1037/NATO/01, 



